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The short time behavior of a disturbed system is influenced by off-shell motion 
and best characterized by the reduced density matrix possessing high energetic 
tails. We present analytically the formation of correlations due to collisions in an 
interacting Fermionic system with and without initial correlation. After this short 
time regime the time evolution is controlled by small gradients. This leads to a 
nonlocal Boltzmann equation for the quasiparticle distribution and a functional 
relating the latter one to the reduced density matrix. The nonlocalities are pre- 
sented as time and space shifts arising from gradient expansion and are leading to 
virial corrections in the thermodynamical limit. 



1 Short time regime 

The generalization of Boltzmann equation towards dense interacting systems 
is a still demanding and unsolved task. A huge variety of different attempts 
can be found in literature to incorporate modificaiiDniLJKhich lead to virial 
corrections in the equation of state, see citations intJfi'tj'QiJ'El. 

These kinetic equations describe different relaxation stages. During the 
very fast first stage, correlations imposed by the initial preparation of the sys- 
tem are decavingLl. During this stage of relaxation the quasiparticle picture 
is established l3'El. After this very fast process the second state develops dur- 
ing which the one-particle distribution relaxes towards the equilibrium value 
with a relaxation time. During this relaxation state the virial corrections are 
established and can be consistently described by a nonlocal Boltzmann kinetic 
equationB. We will present results for both stages here. 

The formation of correlations is connected with an increase of the kinetic 
energy or equivalently the build up of correlation energy. This is due to rear- 
rangement processes which let decay higher order correlation functions until 
only the one - particle distribution function relaxes. Because the correlation 
energy is a two - particle observable we expect that the relaxation of higher 
order correlations can be observed best within this quantity. Since the total 
energy of the system is conserved the kinetic energy increases on cost of the 
correlation energy i?corr(i) and will be calculated from kinetic equations. 

The kinetic equation for the reduced density distribution or Wigner func- 
tion in Born approximation including memory effects but no damping is called 
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Levinson equation and reads in spatial homogeneous media B't3E 



[k-qf ip + q)^ 



{p{T,k~q)p{T,p + q) 



X (1 - p{t, k) - p{t,p)) - pir, k)p{T,p){l - p{T,p + q) - p{t, k - q))} (1) 

with the spin (isospin,...) degeneracy s. The Wigner distribution function is 
normahzed to the density as s J j^jsf{p) — n. The memory effect is con- 
densed in the exphcit retardation of the distribution function. This would lead 
to gradient contributions to the kinetic equation which can be shown to be re- 
sponsible for the formation of high energetic tails in the distribution function 
E3. This will be established on the second stage of relaxation and will lead to 
virial corrections in chapter || 

The second effect is contained in the energy broadening or off-shell behav- 
ior in ([^). This is exclusively related to the spectral properties of the two- 
particle propagator and therefore determined by the relaxation of two-particle 
correlation. On this time scale the memory in the distribution functions can be 
neglected but we will keep the spectral relaxation in the off-shell cos-function 
of The resulting expression for (|l]) describes then how two particles cor- 
relate their motion to avoid the strong interaction regions. This very fast 
formation of the off-shell contribution to Wignpr's distribution has been found 
in numerical treatments of Green's functions Eja. 

Starting with a sudden switching approximation we consider Coulomb in- 
teraction and during the first transient time period the screening is fopaed. 
This can be described by the non-Markovian Lenard - Balescu equationll3 in- 
stead of the static screened equation (|l|) leading to the dynamical expression 
of the correlation energy. To demonstrate its results and limitations, we use 
Maxwellian initial distributions with temperature T neglecting degeneracy. 

From (|l|) 
plasma 



we find with gj-Ecorr = - Z]q / 



dk 
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where we used z 
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t^ 



Up\l L Xt-^j! and Z\ ^l/ ^^rp. 

alytical quantum result of the formation of correlation for statically as well 
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This is the an- 
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as dynamically screened potentials. jlEhc long time limit of (|^) leads to the 
quantum MontroU correlation energyE3 for the dynamical case. 

The situation of sudden switching of interaction may be considered as an 
artificial one. In the simulation experiment we have initial correlations which 
are due to the set up within quasiperiodic boundary condition and Ewald 
summations. Therefore we^XjOusider now initial correlations which lead to 
an additional collision termEj-CJO. This term can be written into the form 
of with the only difference that the occuring Wigner functions are not 
time dependfint but the initial ones and the interaction Vq represents initial 
correlation£Il. The additional collision term, Iq, cancels therefore exactly the 
Levinson collision term ([^) in the case that we have initially the same inter- 
action as during the dynamical evolution (Vq = V) and if the system starts 
from the equilibrium p{t) = pQ. We choose a model interaction of Debye po- 
tential Vi{q) ~ Ame^ l\q^ -|- k|] with fixed parameter Ki = Kjy and for the initial 
correlations Ki = kq to obtain the kinetic energy from the kinetic equations 

Ekin{t) — -Etotal — Einit{t) — Econ{t). (3) 

For the classical limit the time dependent coUisional energy reads 



nT Ax 



-dy{yT{y))y=^r, (4) 



where J-{y) = 1 — e2^^erfc(y), t — tujp/\/2, x = hd/k and = Ane^n/T = 
LUpT/m. The plasma parameter is given as usually by F = ^^j^' '^^^''^ ^.e = 
(jfn)^^'^' ^^S- 0jr-4PPer panel, we compare the analytical results of (P 
with MD simulations t3 using the Debye potential Vi as bare interaction. The 
evolution of kinetic energy is shown for three different ratios x. The agreement 
between theory and simulations is quite satisfactory, in particular, the short 
time behavior for x = 2. The stronger initial increase of kinetic energy observed 
in the simulations at x = 1 may be due to the finite size of the simulation box 
which could more and more influence the results for increasing range of the 
interaction. Now we include the initial correlations which leads to 

iginit(^) \/3F3/2 

^ ^ ^^[xJ'ixT) - XoT{xqt)\, 5 

ni 2[Xq — X'') 

where xq = kq/k characterizing the strength of the initial correlations with 
the Debye potential Vq which contains kq instead of kjj. Besides the kinetic 
energy ^ from initial correlations, the total energy -Etotai now includes 
the correlation energy at the beginning which is the total energy content of 
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Figure 1: The formation of correla- 
tion energy -Ecorr = -Etotal - £^init - 
^coll = ^kin in a plasma with Debye 
interaction Vi. The upper paiicl com- 
pares the analytical pesults (U) with 
MD simulations fromtfl for three dif- 
ferent ratios of kjj to the inverse Debye 
length X = kd/k. In the lower panel 
we compare theoretical predictions for 
the inclusion of Debye initial corre- 
lations characterized by xq = kq/k 
where x = kjj / k = 1. 



the system and can be calculated from the long time limit of (|j) leading to 

i^total V3r3/2 



nT 2{x + xa)' 



(6) 



The result (||) is seen in Fig. ||, lower panel. We observe that if the initial 
correlation is characterized by a potential range larger than the Debye screen- 
ing length, xq < 1, the initial state is over-correlated, and the correlation 
energy starts at a higher absolute value than without initial correlations re- 
laxing towards the correct equilibrium value. If, instead, xq = 1 no change of 
correlation energy is observed, as expected. Similar trends have been observed 
in numerical solutionsll3. 

The characteristic time of formation of correlations at high temperature 
limit is given by the inverse plasma frequency Tc ^ — — This confirms 

the numerical finding ofQ that the correlation or memory time is proportional 
to the range of interaction. In the low temperature region, i.e., in a highly 
degenerated system 2> T, one finds that the build up time is the inverse 
Fermi energy, h/ef in agreement with the quasiparticle formation time known 
as Landau's criterion. Indeed, the quasiparticle formation and the build up of 
correlations are two alternative views of the same phenomenon. 



2 Virial regime 

Now we proceed and investigate how the relaxation into equilibrium is per- 
formed by the system. We will convince ourselves firstly that the direct use of 
reduced density matrix as done so far for short time regimes is not possible. 
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On larger time scales we expect smoothed time-space gradients and expand 
the Levinson equation for the reduced density matrix up to first order memory 
to obtained from (Q) 

— p = Isoltz + -Q^lcoiv (7) 

where besides the usual Boltzmann term, IboUzj there appears an additional 
ofF-shell term, Xcom which describes the correlation in the reduced density 
matrix. For the density balance we see from (j^) that the correlation part from 
the righL.idjdp_/ dpTcorr — i^c which is known as the Beth-Uhlenbeck equation 
of state Ed'Eirca cancels exactly with the left hand side of according to the 
definition J dpp = Uf + ric- So we have to face the exact cancelation of 
in the balance equation and the incorrect conservation of n/. We will show 
that the virial corrections and the correlated density will appear instead from 
internal gradients. That the correlated part Xcorr does not provide the virial 
corrections but a relation between reduced density matrix and quasiparticle 
distribution has not been noted in the literature. 

The second puzzle is a more obvious hint that the Levinson equation fails 
to describe the long time evolution of a system. We can iterate the time 
derivative inside the collision integral on the right hand side of (0) to create 
an infinite series. Neglecting backscattering terms we estimate 

^/>C[/]W.[/] (8) 
with C, = J2 ^^'^ V is similar to the backfiow (second order response) 

n=0 

and is a positive number but not every time smaller than one. Therefore C 
appears to be divergent at certain momentum situations. Therefore we con- 
clude that the gradient expanded Levinson equation leads not to a stationary 
solution in any case. This has also been observed in—numerical solutions of 
these non-Markovian kinetic equations, see figure 4 ofE3. We will show that a 
proper balance between off-shell terms in the reduced density and the kinetic 
equation for the pole part is necessary to ensure stationary solutions. 

The approach presented now is based on the real-time Greenfunction tech- 
nique. We consider the two independent correlation functions for Fermionic 
creation operators G>{1,2) =< a(l)a+(2) > and G<(1,2) a+(2)a(l) >, 
where cumulative variables means time, space, spin, etc 1 — t,x,s.... The 
time diagonal part of yields the reduced density matrix according to 
p(xi,X2,t) = G^(l, 2)(j=t2. For the-jlatter one the time diagonal Kadanoff 
and Baym equation of motion readsc3 
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t 

i[GQ^ , p]{xi, X2,t) — J dt' j dx \p'^ {xi,x ,t^t')T?' {x ,X2Tt' ,t) 

to 

x', t, t')G^ {x' , X2,t' ,t) — {xi,x' , t, i')S^ {x' , X2,t' , t) 
-S<(xi, x', t, t')G> {x', X2, t', t)] (9) 

with the Hartree- Fock drift term Gq^. The right hand side contains a non- 
Markovian colhsion integral and can be considered as a precursor of Levinson 
equation (|^). Using gradient approximation we obtain a kinetic equation of 
the structure and see that the time derivative of the Wigner function on 
the left side of combines with the time derivative of the off-shell part of 
right hand side into a time derivative of a quasiparticle distribution / 

p../-/^-^.<M (10) 

J ZTT Ul ~ € 

with V' the derivative of the principal value. The wave function renormaliza- 
tion has been abbreviated as z = (\ — d^a^^^ k, l+9(^tT which can be confirmed 
using the spectral decomposition <7(lo) — J ^^ ^^^, 7(0;') with the imaginary 
part 7 = (T^ + (T^. We have fulfilled the task and give with ([l^) a connection 
between reduced density matrix p and the quasiparticle distribution /. We 
extrapolate the ansatz also for the correlation functions 

^<f,.,\ _ n^^Xf ^\ I 1 ^ / \ ^^^<i 



g<{Lu) = 2nzS{u;-e)i^ ^ j - -—a< (u) . (11) 

The spectral identity a = + proofs that this ansatz is consistent with 
the extended quasipa rticle pipti irn w h i ch is obtained for small imaginary parts 



of the self-energy 7By ElHamfe 



2tt6(lu - e(prt)) , , d V 

The quasiparticle energies e{prt) are a solution of the dispersion relation uj — 
^ — a{mdrt) =0. It is noteworthy to remark that (|lj) fulfills the spectral 
sum ruleEZI and the energy weighted sura ruleO. The limit of small scattering 
rates has been first introduced by Craiga. An inverse functional f[p] has been 
constructed in cil. . Fpi: equilibrium nonideal plasmas this approximation has 
been employed byEitcZI and under the name of the generalized Beth-Uhlenbeck 
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approach has been used_,by u in nuclear matter for studies of the correlated 
density. The authors inE3 have used this approximation with the name extended 
quasiparticle approximation for the study of the mean removal energy and high- 
momenta tails of Wigner's distribution. The non-equililirium form has been 
derived finally as modified Kadanoff and Baym ansatzi3. Using this ansatz 
( prj ) from the Kadanoff-Baym equations (||) the known Landau-Silin equation 
for the quasiparticle distribution / follows 

^/ + VfeeV,/ - V.eVfe/ = z {{1 - f)a< - fa>) . (13) 



We repeat that the Landau-Silin equation (^^ is coupled with a functional 
that specifies a relation between the quasiparticle distribution / and Wigner's 



function p via ( 10 ) . Using different approximations for the self energy we obtain 
all known kinds of kinetic equations with the generalization that the internal 
gradients of collision integrals will yield the nonlocal or virial corrections. 

The quality of the extended quasiparticle approximation ( p^ ) can be seen 
in the figure of the contribution by P. Lipavsky. One can see that the off-shell 
contribution given by the difference between the Wigner and the Fermi-Dirac 
distributions is not small, in particular at the high momenta region where 
the power-law off-shell tails always dominate over the exponentially falling 
quasiparticle distribution. Formula (|2|) provides inevitable and sufficiently 
precise off-shell corrections. 

In order to show the nonlocal character of the kinetic equation ( p^ ) explic- 
itly, we choose as an approximation for the self energy the ladder or T-matrix 
approximation appropriate for dense interacting systems of short range poten- 
tials 

I]< (1, 2) = T^(l, 3; 5, 6)r^(7, 8; 2, 4)G> (4, 3)G< (5, 7)G< (6, 8). (14) 

For non-degenerate systems, the gradient expansion has been carried through 
by Barwinkela. One can see in Barwinkel's papers, that the scattering integral 
is troubled by a large set of gradient corrections. This formal complexity seems 
to bnJjjie main reason why most authors either neglect gradient coxmctions at 
all c3cB or provide them buried in multi-dimensional integrals eJ'l3. For a 
degenerate system, the set of gradient corrections to the scattering integral is 
even larger than for rare gases studied by Barwinkel, seeB. 

The quasiclassical limit with all linear gradients kept is a tedious but 
straightforward algebraic exercise. It results in one nongradient and nineteen 
gradient terms that are analogous to those found within the chemical physics 
All these terms can be recollected into a nonlocal and noninstantaneous 
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scattering integral ttiat has an intuitively appealing structure of a nonlocal 
Boltzmann equation □ 



dfi ^ dei dfi dei dfi ^ 
dt dk dr dr dk ^ 



dpdq 
(27r)5 



S (ei + S2 - es - Si + 2Ae) Z1Z2Z3Z4, 



%b[ ei+e2-AE,k--^,p-^,q,r-Ar,t-^ 



/3/4(l - /l) (1 - /2) - (1 - h) (1 - /4)/l/2 



\_ldA2 de2dA2^ 




with Enskog-type shifts of argumentsB: /i = fa{k,r,t), /a = r- A2, i), 
/3 = faik~q-A^r-A3,t-At), and U = fbip + q- Ak ,r~ A^^t- At). In 
agreement withLJa, all gradient corrections result proportional to derivatives 
of the scattering phase shift = Im ln'7^^(ri, k,p, q, t, r), 



A, = 



A, = 



d(j) d(j) 



dp dq dk 



A3=-^ 

^ dk 



£3+^4 



^ dk ^ dq 



£3+54 



dn 



A, 



£3 +£4 



190 

2 at 



£3 +£4 



190 
2 9^ 



£3 +£4 



and Ar 
shell 



(16) 

i(A2 + A3 + A4). After derivatives, A's are evaluated at the energy 
£4- 



4 

£3 



The A's in arguments of distribution functions in (15) remind non-instant 
and non-local corrections in the scattering-in integral for classical particles. 
The displacements of the asymptotic states are given by A2,3,4. The time 
delay enters in an equal way the asymptotic states 3 and 4. The momentum 
gain Ak also appears-only in states 3 and 4. Finally, there is the energy gain 
which is discussed inca. We remind that the non-localities should be viewed 
as a compact form of gradient corrections. 

Despite its complicated form it is possible to solve this kinetic equation 
with standard Boltzmann numerical codes and to implement the shifts l3. 
Themfore we have calculated the shifts for different realistic nuclear poten- 
tials EZI. The numerical solution of the nonlocal kinetic equation has shown an 
observable effect in the dynamical particle spectra of around 10%. The high 
energetic tails of the spectrum is enhanced due to more violent collisions in 
the early phase of nuclear collision. Therefore the nonlocal corrections lead to 
an enhanced production of preequilibrium high energetic particles. 
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3 Thermodynamic properties 



The meaning of nonlocal shifts can be best seen on thermodynamic observables 
like density of particles a, density of energy £, and the stress tensor Jij 
which conserve within the nonlocal and noninstantaneous kinetic equation. (^5|). 
Integrating ( |l5|) over momentum k with factors £i, k and unity one findsQ that 
each observable has the standard quasiparticle part following from the drift 




and the A-contribution following from the scattering integral 

a,b ^ ' 

a.b 

A„.^/^Pi„ (18) 

where P = \Tab\''2Tr5{ei+e2-S3~ei)hf2il- h- U). 

The density of energy £ — E'^^ + A£ alternatively results from Kadanoff 

and Baym formula, £ = / ^ / + ^) G<(w, fc, r, t), with G< in 

the extended quasiparticle approximation (pi]). Its complicated form, how- 
ever, shows that £ cannot be easily inferred from an_evcntual experimental fit 
of the kinetic equation as it has been attempted intS. The con^jCrvation of £ 
generalizes the result of Bornath, Kremp, Kraeft and Schlangescl restricted to 
non-degenerated systems. The particle density Ua = n'^^ + Ana is also obtained 
from (|l^) via the definition, Ua = J '^j0ysG<. This confir ms that the ex- 
tended quasiparticle approximation is thermodynamically consistent with the 
nonlocal and noninstantaneous corrections to the scattering integral. 

For equilibrium distributions, formulas (^|) and provide equations of 
state. Two known cases are worth to compare. First, the particle density rL^ff 
nJJP -|- Aua is identical to the quantum Beth-Uhlenbeck equation of state EatI, 



9 



where nJJP is called the free density and Aua the correlated density. Second, the 
virial correction to the stress tensor has a form o£ the collision flux contribution 
known in the theory of moderately dense gasesu'B. 

4 Summary 

While the short time behavior is described by off-shell transport condensed 
in the reduced density matrix, the later stage of evolution is described by a 
nonlocal Boltzmann equation. This equation is derived in the quasiclassical 
limit and leads to consistent thermodynamics in equilibrium including binary 
correlations. The latter ones represent the second virial coefficients in the low 
density limit. The presented kinetic equation unifies achievements of dense 
gases and the kinetic theory of quasiparticle transport in quantum systems. 
The equation has been shown applicable in recent simulation codes. 
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